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Abstract 

The completeness of the group classification of systems of two linear second-order 
ordinary differential equations with constant coefficients is delineated in the paper. 
The new cases extend what has been done in the literature. These cases correspond 
to the type of equations where the commutative property of the coefficient matri- 
ces with respect to the dependent variables and the first-order derivatives in the 
considered system does not hold. A discussion of the results as well as a note on 
the extension to linear systems of second-order ordinary differential equations with 
more than two equations are given. 
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1 Introduction 



In this paper we consider the complete group classification of systems of two 
linear second-order ordinary differential equations with constant coefficients. 
Systems of second-order ordinary differential equations arise in the modeling 
of physical phenomena in the areas of physics, chemistry and mathematics. 
Of interest to the study of such systems is their symmetry properties. The 
founder of symmetry analysis of scalar ordinary differential equations, Sophus 
Lie [T], gave a complete group classification of a scalar ordinary differential 
equation of the form 

y" = f(x,y). 

L.V.Ovsiannikov j2j later performed Lie's classification in a different way. This 
classification was obtained by directly solving the determining equations and 
using the equivalence transformations. The same approach was applied in 
[3] for the group classification of more general types of equations of the form 
y" = P${x, y; y'), where Ps{x, y; y') is a polynomial of third degree with respect 
to the first-order derivative y' . In the general case of a scalar ordinary differ- 
ential equation y" = f(x,y,y') the application of the method that involves 
directly solving the determining equations gives rise to some difficulties. The 
group classification of such equations [I] is based on the enumeration of all 
possible Lie algebras of operators acting on the plane (x,y). Lie [5] gave the 
classification of all dissimilar Lie algebras (under complex change of variables) 
in two complex variables. Later on, in 1992, Gonzalez-Lopez et al. ordered the 
Lie classification of realizations of complex Lie algebras and extended it to 
the real case [6]. A large amount of results on the dimension and structure of 
symmetry algebras of linearizable ordinary differential equations is well-known 

(see |H7i5ignniini2] ). 

There are several papers where the group classification of systems of second- 
order ordinary differential equations was considered. The literature has dealt 
extensively with symmetry properties of a scalar ordinary differential equation 
but the group classification of a system of even two linear second-order equa- 
tions with constant coefficients is not complete. In recent works [T2|13|14|15|16] 
the authors focused on the study of systems of second-order ordinary differ- 
ential equations with constant coefficients of the form 

y" = My, (1) 

where M is a matrix with constant entries. In the present paper it is shown 
that these types of systems do not exhaust a set of all systems of linear second- 
order ordinary differential equations with constant coefficients. 

In the general case of systems of two linear second-order ordinary differential 
equations, the more advanced results are obtained in [T7], where using the 
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canonical form, 

[ a(x) b(x) \ 

y = y, 

\ c(x) — a(x) J 

presented several representatives of nonequivalent classes. However the list of 
all distinguished representatives of systems of two linear second-order ordi- 
nary differential equations was not obtained in [17]. Recently the complete 
group classification of linear systems of two second-order ordinary differential 
equations has been studied in [IB] . 

Here the complete classification of linear systems of two second-order ordinary 
differential equations with constant coefficients is presented. 

As far as we are aware the results found here are new and have not been 
reported in the literature. 

The paper is organized as follows: In the first part we present the simplification 
of a linear system of second-order equations with constant coefficients. This is 
then followed by a preliminary study of linear systems of second-order ordinary 
differential equations with constant coefficients. The later part is followed by 
a complete treatment of linear systems of second-order ordinary differential 
equations with constant coefficients as well as the discussion of the results and 
conclusion in the final part of the paper. 



2 Simplification of a system of linear equations 

Let a linear system of second-order equations with constant coefficients be 
given by: 

y " = Ay' + By + f, (2) 

where A and B are constant matrices. Using a particular solution y p , one can 
reduce (j2J) to the homogeneous system 

y" = Ay' + By. (3) 
2. 1 Linear change of the dependent variables 
Applying the change, 

y = Cx, 

where C = C(t) is a nonsingular matrix, system ([3]) becomes 

x" = Ax! + fix, (4) 



3 



where 

A = C-\AC - 2C') , B = C~ X {BC + AC' - C") . 
If one chooses the matrix C(t) such that 

C = l -AC, 

then 

B = C~ l (BC + l -A 2 C - ^A 2 c\ = C- 1 (b + ^A 2 ^j C. 

The existence of the nonsingular matrix C(t) is guaranteed by the existence 
of the solution of the Cauchy problem 

C'= l ~AC, C(0) = E, 

where E is the unit matrix. This solution defines the matrix 

C(t) = eH 
which commutes with the matrix A: 

AC = CA. 

Let us study the case when the matrix B is constant. First of all notice that 

A (c- 1 ) = -c- l cc-\ 

dt v / 

Then one has 

i* = i { c ~ l i B + c ) 

= -C- l ACC- 1 (B + \A 2 ) C + C- 1 (B + \A 2 ) AC 
= -C- l A (B + \A 2 ) C + C- 1 (B + \A 2 ) AC = C~ l {BA - AB)C. 

Hence, the matrix B is constant if and only if the matrices A and B commute 

AB = BA. (5) 

The papers cited in the literature [T2fl3|ll4|ll5|ll6] . where the group classifica- 
tion of systems of linear equations with constant coefficients was considered, 
studied systems with constant matrix B. Thus, even for systems of linear 
equations with constant coefficients there is no complete study of the group 
classification. The present paper fills this niche. 
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2.2 Example 



Let us choose the matrices 



then 

Notice that 
Then 



Hence, 

where <p(2t) 
that 





foi] 




foi\ 


A = 




, B = 






{ 10) 




[oo) 



AB-BA 



-1 
1 



n sA 



A 2 = E. 

E + sA+^A 2 + ^A :i + ^A 4 + ... 

(£. _L §1 _L I ^ 2 " +1 "\ P 

" V2! ^ 4! ^ ••• ~r (2n+l)!"V 

+ (^+^ 3 + - + (SrV^ +1 + -) 

= I ((e s + e" s ) £ + (e s - e- s )) A. 



C(t) = <p(t)E + if>(t)A, 
\ (e* + e~') and ^(2*) = § (e* — e - *). One can check directly 



i = Q £ = 1 f (1 + ^(2*)) (2 + ^(2t)) 
4^(2-^(2*)) (1-^(20) 
It is obvious that the matrix B is not constant. 



3 Preliminary study of linear systems of second-order ordinary 
differential equations with constant coefficients 



Let us consider a linear system of second-order ordinary differential equations, 

y" = V + £y, (6) 

where A and B are matrices with constant entries. For commuting matrices 
AB = BA this system is reduced to the system of the form 

y" = My, (7) 
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which is completely studied in [T2"fl3|ll4|ll5|ll6] . The case of noncommutative 
matrices is considered here. 

3.1 Equivalence transformations 

The general solution of the determining equations of the equivalence Lie group 
defines the equivalence Lie group, which corresponds to the generators 

Xl = d x , X e 2 = yd y + zd z , 

X| = zd z — a\ 2 d ai2 + a2id a21 — b\ 2 d bl2 + b2idb 21 , 
X% = yd z - a 12 d ail - (a 22 - a u )d a21 + a 12 d a22 

-h 2 d bll - (622 - b u )d b21 + b 12 d b22 , 
X§ = zd y + a 2 id an + (a 22 — an)d ai2 — a 2 \d a22 

+b 2 id bll + (b 22 - b n )d bl2 - b 21 d b22 , 
Xg = xd x — aud au — a 12 d ai2 — a 2 \d a21 — a 22 d a22 

-2b u d bll - 2b 12 d bl2 - 2b 21 d b , 21 - 2b 22 d b22 , 
X 7 e = x{yd y + zd z ) + 2d ail + 2d a22 

—and bu — ai 2 d bl2 — a 2 \d b21 — a 22 d b22 . 

The transformations related with the generators X|, X|, X% and X§ corre- 
spond to the linear change of the dependent variables y = Py with a constant 
nonsingular matrix P. The generators Xf and Xg define shifting and scaling 
x. The transformations corresponding to the generator Xj define the change 

y = ye TX , z = ze TX , a u = a u + 2r, a 12 = a 12 , 
a 2 \ = 021) o 2 2 = a 22 + 2r, bn = bn — «iiT — r 2 , 
612 = fen - a 12 T, b 21 = b 21 - a 21 r, b 22 = b 22 - a 22 r - r 2 , 

where r is a group parameter. In matrix form the last set of transformations 
is rewritten as 

y = e TX y, A = A + 2rJ, B = B-tA- t 2 I, (8) 
where I is the identical matrix. 



6 



3.2 Simplification of the matrix A 



Let us apply the change y = Py, where P is a nonsingular matrix with 
constant entries 

/, 

P = 

\P21 P22 



Pn P12 



Ay 1 + By, 



Equations (EJ) become 

y 

where 

A = PAP-\ B = PBP-\ 

This means that the change y = Py reduces equation ([6]) to the same form 
with the matrices A and B changed. Using this change, the matrix A can be 
presented in the Jordan form. For a real-valued 2x2 matrix A, and if the 
matrix P also has real-valued entries, then the Jordan matrix is one of the 
following three types, 



Ji 





a 1 
a 



(9) 



where a, b and c > are real numbers. By virtue of transformations ([H]) one 
can assume that a = 0. Notice also that using the dilation of x, one can reduce 
c to 1. 



In these cases one obtains the following three cases: 
A = J 1[a=0 -AB-BA-- 

A = J 2 \ a =o,c=i : AB — BA = 



A = J 3la=0 :AB-BA 



where 6^ are entries of the matrix B. 



-b 12 b 
b 21 b 

O12 + b 21 ) {b 22 - b u ) 
(b 22 - b n ) -(612 + 621, 
hi (622 - b n ) 
-b 21 



This then gives the following for noncommutative matrices: 

(a) in the case A = J x one can assume that b\ 2 b ^ 0; 

(b) in the case A = J 2 one has to assume that (b u + b 21 ) 2 + {b 22 — bn) 2 ^ 0; 

(c) in the case A = J 3 one has to assume that b 2 2 + ib 22 — bn) 2 ^ 0. 
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3.3 Case A = Ji 



In this case the matrix A is given as follows: 

A = 




Let us denote b = 4A, where A ^ 0. The admitted generator has the form 

X = C\X\ + C 2 X 2 , 
and the determining equations are reduced to the study of the equations 

CiCftfi - 2b 2 u b 22 + 76?! A 2 + 2&U&U&21 + &ii&la " A 2 
-566uA 4 - 2b 12 b 2l b 22 - b 2 22 X 2 + 8& 22 A 4 + 48 A 6 ) = 0, 



(10) 



= 0, (11) 

(24A 4 + 146 22 A 2 - 6&uA 2 - 2b u b 22 + b 12 b 21 + 26 2 2 )d = 0, (12) 

C 2 b 2l = 0, (13) 

(Ab 22 + 15A 2 )C 2 = 0, (4& n - A 2 )C 2 = 0, (14) 



where 



hi = b u + & 22 + 2A 2 , h 2 = b 22 - b u + 4A 5 
X x = e- 2Xx (h 2 (d x - \{yd y - zd z )) - 2\b 12 zd, 



V! 1 



X 2 = e~ Xx (2d x - X(yd y - 3zd z )) . 



Let us consider equations which admit the generator X 2 . In this case C 2 ^ 
and from equations (1T31) - (1T4I) one finds that 

6 n = A 2 /4, 6 22 = -15A 2 /4, b 21 = 0, (15) 

or the matrix 

1 ( A 2 46 12 \ 

B = -A • ( ie ) 

4 ^ -15A 2 ) 

Notice that in this case h 2 = 0. Since in this case equations (TTU ]) - (fr2|) are also 
satisfied, then there are two admitted generators 

X x = e- 2Xx zd y , X 2 = e~ Xx (2d x - \{yd y - 3zd z )) . (17) 



If the matrix B is not of the form (1161) . then C 2 = 0, and for the existence of a 
nontrivial admitted generators one needs to require that C\ 7^ 0. Considering 
equation (fTTj) . first assume that h\ = 0. Then equations ffTUj) and (fl2"]l give 

611 = A 2 , 6 22 = -3A 2 , 6 2 i = 0. 
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These conditions provide that h 2 = 0. Thus the case Y12 = is the general 
case. 

Assume now that h 2 = 0. The general solution of equations ffT0]) - ffT2]) is 

b u = b 22 + 4\ 2 , 6 2 i = 0, 

or the matrix 

_ I 622 + 4A 2 612 
V 622 
In this case equations ffl4|) are reduced to the equation 

C 2 (4622 + 15A 2 ) = 0. 

Thus, if 46 2 2 + 15 A 2 = 0, then one comes to the studied case (1T51) . where 
the admitted generators are f)17p . If 46 2 2 + 15A 2 7^ 0, then there is the only 
admitted generator X\. 

In the summary the result is the following. 

There are nontrivial admitted generators only if the matrix B has the form: 



B 



622 + 4A 2 612 
60 



\ u o 22 

If 4622 + 15A 2 7^ 0, then there is the only admitted generator 

X x = e~ 2Xx zd y . 

The extension of the Lie algebra, defined by the generator X\, only occurs for 
622 = — 15A 2 /4. In this case the equations admit two generators 

X x = e- 2Xx zd y , X 2 = e~ Xx (2d x - X(yd y - 3zd z )) . 
3.4 Cases A = J 2 and A = J 3 



In these cases 

/ 

A = 

V 



-10/ loo 



respectively. Calculations show that there are no extensions of the trivial gen- 
erators in these cases. 
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4 Conclusion 



This paper completes the study of the group classification of linear systems 
of two second-order ordinary differential equations with constant coefficients. 
In the literature the studied case reduced to y" = My where M is a ma- 
trix with constant entries. This study can be extended to linear systems of 
second-order ordinary differential equations with more than two equations. 
For noncommutative matrices we found new cases where there are additional 
generators in addition to the generic ones, d x and yd y + zd z . This implies that 
when considering the classification of such systems the first step should be to 
check whether the matrices A and B commute in (jH]). If they commute then 
the case is equivalent to the case ([1]) studied in the literature where M is a 
constant matrix whilst if A and B do not commute then the study in this 
paper applies and can be extended to deal with linear systems of second-order 
ordinary differential equations with more than two equations. 
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